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The origins of the anomalous temperature dependence of magnetocrystalline anisotropy in 
(Fei_ccCox)2B alloys are elucidated using first-principles calculations within the disordered local 
moment model. Excellent agreement with experimental data is obtained. The anomalies are asso¬ 
ciated with the changes in band occupations due to Stoner-like band shifts and with the selective 
suppression of spin-orbit “hot spots” by thermal spin fluctuations. Under certain conditions, the 
anisotropy can increase, rather than decrease, with decreasing magnetization due to these peculiar 
electronic mechanisms, which contrast starkly with those assumed in existing models. 


Magnetocrystalline anisotropy (MCA) is one of the key 
properties of a magnetic material [1]. Understanding of 
its temperature dependence is a challenging theoretical 
problem with implications for the design of better mate¬ 
rials for permanent magnets [ 2 ], heat-assisted magnetic 
recording [3], and other applications. While the MCA 
energy K usually declines monotonically with increasing 
temperature as predicted by simple models [4], in some 
magnets it behaves very differently and can even increase 
with temperature. Such anomalous K{T) dependence 
makes some materials useful as permanent magnets and 
can potentially facilitate specialized applications. 

Well-known anomalies in the temperature dependence 
of MCA include spin reorientation transitions (SRT) in 
cobalt [5] and MnBi [ 6 ], which have been attributed to 
thermal expansion; an SRT in gadolinium, which may be 
due to higher-order terms in MCA [7]; SRT in R 2 Fei 4 B 
hard magnets [ 8 ] due to the ordering of the rare-earth 
spins at low T; and SRT in thin films [9, 10] associated 
with the competition between the bulk and surface con¬ 
tributions to MCA. Competition between single-site and 
two-site MCA can also lead to an SRT [ 11 ]. 

MCA in metallic magnets is rarely dominated by the 
single-ion mechanism leading to the K oc dependence 
on the magnetization [4]. For example, two-ion terms in 
3d-bd alloys like FePt modify this dependence to iF oc 
[ 12 , 13]. Clear understanding of the anomalous 
temperature dependence of MCA has been so far limited 
to the cases when competing contributions to MCA can 
be sorted out in real space, such as, for example, bulk and 
surface terms in thin films. In contrast, understanding of 
MCA in itinerant magnets usually requires a reciprocal 
space analysis [14]. 

One such system is the disordered substitutional 
(Fei_a^Coa^) 2 B alloy, which exhibits three concentration- 
driven SRTs at T = 0, a high-temperature SRT at the 
Fe-rich end, and a strongly non-monotonic temperature 
dependence at the Co-rich end with a low-temperature 
SRT [15, 16]. The SRT’s at T = 0 were traced down 
to the variation of the band filling with concentration 


combined with spin-orbital selection rules [16]. Here we 
elucidate the unconventional mechanisms leading to the 
spectacular anomalies in the temperature dependence of 
MCA in this system and show that they stem from the 
changes in the electronic structure induced by spin fluc¬ 
tuations. We will see that under certain conditions MCA 
can increase, rather than decrease, with decreasing mag¬ 
netization due to these mechanisms. 

Our calculations employ the Green’s function-based 
linear muffin-tin orbital method [17] with spin-orbit cou¬ 
pling (SOC) included as a perturbation to the potential 
parameters [16, 18]. Thermal spin fluctuations are in¬ 
cluded within the disordered local moment (DLM) model 
[19, 20], which treats them within the coherent potential 
approximation (CPA) on the same footing with chemical 
disorder. The DLM method has been previously used 
to calculate the K(T) dependence in systems like FePt 
[ 21 , 22 ] and YC 05 [23]. Although K(T) in these metals 
does not follow the Callen-Callen model [4] designed for 
materials with single-ion MCA, it still decreases mono¬ 
tonically. In contrast, we will see that the changes in 
the electronic structure with temperature lead to strong 
anomalies in (Fei_a^Coa^) 2 B. Our implementation of the 
DLM method is described in Ref. 24. (See Supplemental 
Material [25] for additional details.) 

Apart from the inclusion of spin disorder, the com¬ 
putational details are similar to Ref. 16. In particular, 
the large overestimation of the magnetization in density- 
functional calculations for C 02 B ( 1.1 fiB compared to 
experimental 0.76 i^b per Co atom) is corrected by scal¬ 
ing the local part of the exchange-correlation field for Co 
atoms by a factor 0.8 at all concentrations. This treat¬ 
ment is consistent with spin-fluctuation theories showing 
that spin fluctuations tend to reduce the effective Stoner 
parameter [26, 27] and allows us to take into account the 
resulting changes in the electronic structure. 

Magnetism in (Fei_a^Coa^) 2 B alloys is much more itin¬ 
erant compared to systems like FePt; the spin moments 
of Fe and, especially, Co atoms are not rigid in density- 
functional calculations. To implement spin disorder 
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within the DLM method, we make a simple assumption 
that the spin moments of both Fe and Co at finite T can 
be taken from the ferromagnetic state at T = 0. This as¬ 
sumption is based on the expectation that thermal spin 
fluctuations to a large extent restore the “soft” spin mo¬ 
ments [26]. On the other hand, the variation of the elec¬ 
tronic structure with T should not be very sensitive to 
the details of the spin fluctuation model. For simplicity, 
a similar approach is used for the (Coi_a,Nia,) 2 B system, 
including the small spin moments on the Ni atoms. 

The distribution functions for spin orientations are 
taken in the Weiss form: Py{0) oc exp(aiy cos6>), where 
0 is the angle made by the spin with the magnetization 
axis, and v labels the alloy component. The temperature 
dependence of the coefficients is determined using the 
calculated effective exchange parameters, as explained in 
the Supplemental Material [25]. Fermi-Dirac smearing 
is neglected, because the effects of spin fluctuations are 
overwhelmingly stronger. 

The results of K{x^T) calculations shown in Fig. 1, 
which were obtained with temperature-independent lat¬ 
tice parameters, are in excellent agreement with experi¬ 
mental data [15]. Both the thermal SRT at the Fe-rich 
end and the non-monotonic temperature dependence at 
the Co-rich end in (Fei_a,Coa,) 2 B alloys are captured (see 
Supplemental Material [25] for a direct comparison). For 
(Coo. 9 Nio.i) 2 B the MCA energy at T = 0 is large and 
negative in agreement with experiment [15], although the 
initial decline of K(T) similar to C 02 B is not observed 
in experiment. The finite slope in the K(T) curves at 
zero temperature is due to the classical treatment of spin 
fluctuations. We have explicitly verified that the effect of 
thermal expansion on K(T) in Fe 2 B and C 02 B is almost 
unnoticeable. 

The effects of spin disorder on the electronic structure 
can be understood from Fig. 2, which shows the partial 
minority-spin Bloch spectral function at x = 0.95 for T = 
0 and T/Tc = 0.7. Here, at the Co-rich end, all bands 
are easily identifiable and relatively weakly broadened at 
TjTc = 0.7. In addition, they are shifted down relative 
to their positions at T = 0, which is a hallmark of an 
itinerant Stoner system. In contrast, at the Fe-rich end 
the bands are strongly broadened by spin fluctuations, 
so that most bands in the 1 eV window below Ep are 
barely visible (see Supplemental Material [25]). The large 
difference in the degree of band broadening between the 
Fe-rich and Co-rich ends is due to the 2.5-fold difference 
in the magnitude of the spin moments. The effect of 
phonon scattering on band broadening in (Fei_a,Coa,) 2 B 
alloys is likely much smaller and is neglected here. 

The usual expectation is that spin disorder should re¬ 
duce MCA as a result of averaging over spin directions. 
Such normal behavior is seen, for example, at x = 0.3 
in Fig. 1. This expectation is violated at many concen¬ 
trations: K{x^T) is non-monotonic with respect to T at 
0 < X < 0.2, 0.5 < X < 0.6, and 0.9 < x < 1; we will 



FIG. 1 . Calculated temperature dependencies of MCA energy 
K in (Fei-ccCox)2B and (Coo.9Nio.i)2B alloys. 

call this behavior anomalous. At x < 0.6 the anomalous 
temperature dependence of A" at a given x follows the 
variation of K with increasing x at T = 0. For exam¬ 
ple, i^(0.2,0) > i^(0.1,0), and i^(0.1,T) anomalously 
increases with T. At x > 0.9 the anomalous variation is 
opposite to the trend in K{x^0) with increasing x. To 
understand this difference, we first need to examine the 
effect of disorder on MCA. 

Fig. 3 compares K{x,0) calculated within the vir¬ 
tual crystal approximation (VCA) with CPA results for 
(Fei_a^Coa^) 2 B [16] and (Coi_a;Nia,) 2 B systems [28]. Note 
that in the (Coi_a;Nia,) 2 B system the spin moments van¬ 
ish near 40% Ni, in agreement with experiment [29]. In 
addition to the MCA energy K, Fig. 3 also shows its 
approximate spin decomposition obtained from the 
SOC energy [16, 25]. Because the 3d shell in this system 
is more than half filled, the variation of MCA with x is 
largely controlled by the term, i.e., by the LzSz mix¬ 
ing of the minority-spin states. Substitutional disorder 
strongly suppresses MCA, an effect that was also found 
in tetragonal Fe-Co alloys [30]. The suppression is due 
to band broadening, which reduces the efficiency of spin- 
orbital selection rules. Importantly, bands broaden at 
different rates; the contributions to MCA from the bands 
that lie close to Ep and broaden strongly are most effec¬ 
tively suppressed. The dispersive majority-spin bands 
are weakly broadened, and hence the term is almost 

unaffected by disorder; in contrast, is strongly re- 
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FIG. 2. Partial minority-spin spectral function for the 
transition-metal site in (Feo.05000.95)26 at (a) T = 0, and 
(b) T/Tc = 0.7. SOC is included, M || z, and energy is in 
eV. Color encodes the orbital character of the states. The 
intensities of the red, blue and green color channels are pro¬ 
portional to the sum of m = ±2 {xy and sum of 

m = dll {xz and yz), and m = 0 {z^) character, respectively. 


duced. We note that although band broadening (and 
thereby MCA) can depend on chemical short-range or¬ 
der, the latter is expected to be negligible in the present 
alloy with chemically similar constituents. 

The strongest suppression of MCA can be expected for 
the “hot spots” appearing when nearly degenerate bands 
dit Ef are split by SOC [14]. A clear example of such 
bands is seen near the F point in Fig. 2a. The effect 
of disorder is further illustrated in Fig. 4 showing the 
spectral function at the F point for two orientations of 
the magnetization at x = 1, 0.9, and 0.8, all at T = 0. 
At X = 1 there is no disorder, and the sharp bands are 
fully split by SOC for M || z. With the addition of Fe, 
the broadening quickly exceeds the original SOC-induced 
splitting, and the effect of SOC is strongly suppressed. 

Disorder has a similar effect on the mixing of electronic 
bands of opposite spin by L^S- and Indeed, 

while in Fig. 2a for T = 0 the anticrossings with the 
majority-spin bands are clearly visible, in Fig. 2b, for 
T/Tc = 0.7, they are almost completely suppressed. 

We now return to the analysis of the anomalous tem- 





FIG. 3. MGA in (Fei_ccGox) 2 B and (Goi_ccNia,) 2 B alloys 
calculated within VGA (empty circles) compared with GPA 
(hlled circles). The spin decomposition is given for VGA. 
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FIG. 4. Spectral functions at me F point at (a) a: = 1, (b) 
X — 0.9, (c) X — 0.8. Solid lines: M || z. Dashed lines: 
M \\ X. A small imaginary part is added to energy to resolve 
the bands in panel (a). 


perature dependence of K. We expect that these anoma¬ 
lies come from the effects of thermal spin fluctuations on 
the electronic structure beyond a simple averaging over 
spin directions. As we saw in Fig. 2, there are two such 
effects in (Fei_a,Coa,) 2 B: reduction of the exchange split¬ 
ting A, and band broadening. The reduction of A shifts 
the minority-spin bands downward relative to just 
as the band filling with increasing x does. Band broad¬ 
ening has a stronger effect on the minority-spin states, 
where Ep lies within the relatively heavy 3d bands, and 
it is particularly important for nearly degenerate bands 
straddling the Fermi level, as we saw in Fig. 4. 

To understand how these effects lead to to the anoma¬ 
lies in K(T)^ it is convenient to examine two quantities. 
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Kf and K^, defined as = J^°{E - Eo)AN^{E)dE, 
where is the Fermi energy in the absence of SOC, and 
ANa- is the difference, between M || x and M || z, in 
the partial density of states for spin a in the global ref¬ 
erence frame. Their sum closely approximates 

K, and their analysis can help identify the contributions 
of different bands to particularly in combination with 
reciprocal-space resolution [16, 25]. 

Fig. 5a shows the temperature dependence of Kcr in 
Fe 2 B. Since the spin-mixing contribution here is 
small (Fig. 3), and provide information similar to 
and at T = 0 while retaining clear meaning at 
finite temperature [25]. We see that decreases quickly 
with increasing T. This happens because the downward 
shift and broadening of the minority-spin bands strongly 
suppress the negative minority-spin contribution to K. 
In contrast, the initial increase in mirrors the upward 
slope of 0 ) as a function of x [16], which occurs as 

the majority-spin bands shift upward relative io Ep with 
decreasing A. At elevated temperatures the majority- 
spin contribution becomes dominant, and K undergoes 
an anomalous sign change, i.e., a spin-reorientation tran¬ 
sition. 




FIG. 5 . Contributions to K in (a) Fe2B and (b) 
(Feo.osCoo.95)26 from different spins and K\). and 
K~ in panel (b): total positive and negative contributions to 
Ka coming from different k points. (Dotted lines show , 

K,-.) 

At the Co-rich end the situation is complicated by the 
presence of large contributions of opposite sign that come 


from the minority-spin states in different regions of the 
Brillouin zone [16]. Near the F point there is a large 
positive contribution from the degenerate bands that are 
mixed by Lz. There is also a large negative contribution 
from the mixing of minority-spin bands of opposite parity 
with respect to reflection, which is distributed over the 
whole Brillouin zone. To help resolve these contributions. 
Fig. 5b for (Fco.osCoo. 95)26 shows, in addition to the 
total positive {K^) and negative {K~) contributions to 
i^cr, which were sorted by wave vector. Fig. 6 displays 
k-resolved on the F MX plane at T = 0 and T/Tc — 
0.7. The bright red ring around the F point in Fig. 6 
is the hot spot coming from the two nearly-denegerate 
bands that are split by SOC (see Fig. 2 a and 4). 

As seen in Fig. 6 , thermal spin disorder strongly sup¬ 
presses the hot spot observed at T = 0: it is strongly 
washed out at TjTc = 0.7, while the contributions from 
other regions decline almost homogeneously. This effect 
is similar to that of chemical disorder (Fig. 4). As a re¬ 
sult, declines faster compared to other contributions 
shown in Fig. 5b, and the negative value of K grows 
anomalously with T. 

Interestingly, while in VCA the maximum in i^(x,0) 
with respect to band filling occurs near x = 0.95 (Fig. 3), 
in CPA there is a cusped maximum exactly in C 02 B. The 
latter is due to the fact that the bands are broadened by 
disorder with any admixture, reducing the positive con¬ 
tribution from the hot spots. This dominant effect of 
disorder explains why, as noted above, the anomalous 
K{T) dependence at x > 0.9 is opposite to the trend 
expected from increasing x, which holds at other con¬ 
centrations. In C 02 B, where the positive contribution is 
at its maximum, both band broadening and decreasing 
A contribute to the anomalous decrease in K(T)^ as the 
nearly degenerate bands broaden and sink below Ep- 


[Ml E 



FIG. 6. Wave vector-resolved (units of meVao, where ao is 
the Bohr radius) on the FMX plane in (Feo.05 Goo. 95)26 alloy 
at T = 0 (upper left) and T/Tc = 0.7 (lower right). 
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In conclusion, we found that the anomalous tempera¬ 
ture dependence of MCA in (Fei_a,Coa,) 2 B alloys is due 
to the changes in the electronic structure induced by spin 
fluctuations. This unconventional mechanism can be har¬ 
nessed in applications where temperature-independent or 
increasing MCA is required. 
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SUPPLEMENTAL MATERIAL 


DESCRIPTION OF SPIN DISORDER 

All calculations were performed using our implementa¬ 
tion [SI] of the disordered local moment (DLM) model for 
partially ordered magnetic states. Thermal spin disorder 
was introduced as follows. Fig. SI shows the effective 
exchange parameters Jg = d‘^EldO‘f^ for both compo¬ 
nents (/i = Fe, Co) calculated at all concentrations using 
the linear response technique [S2]. These parameters de¬ 
scribe the exchange interaction of an atom of a given 
type with the rest of the crystal. Within the mean-field 
approximation (MFA), these data predict the Curie tem¬ 
peratures Tc of pure Fe 2 B and C 02 B to be 1570 K and 
290 K, which can be compared with experimental values 
of 1013 K and 429 K, respectively [S3]. 



FIG. SI. Effective exchange parameters Jq for Fe and Co 
(dashed lines) and calculated Curie temperature Tc (solid 
line) as a function of concentration. Blue triangles: experi¬ 
mental data for Tc [S3]. 

For Fe 2 B the Tc is overestimated by about 35%, which 
can largely be attributed to the neglect of short-range 
order in MFA. Indeed, linear response calculations show 
that the spin of a Fe atom is strongly ferromagnetically 
coupled to seven neighbors, five of which are within the 
layer. Although the couplings to more distant atoms are 
not negligible, they alternate in sign and contribute little 
to the effective exchange parameter. Thus, MFA can be 
expected to overestimate Tc by about as much as it does 
(about 23% [S4]) for the bcc lattice with nearest-neighbor 
exchange (coordination number 8). In contrast, for C 02 B 
the Tc is underestimated by about a factor 1.5. This sit¬ 
uation appears to be similar to the well-known case of 
fee Ni, where Tc calculated from Jo is underestimated 
by nearly a factor of 2 [S5]. The similarity is due to the 
fact that both Ni and C 02 B are strongly itinerant fer- 
romagnets with relatively small local moments. In such 
systems the effects of quantum spin fluctuations become 
significant, which for C 02 B is also reflected in the large 
overestimation of the magnetization in density-functional 
calculations. In addition, the long-wave approximation 
inherent in the linear-response calculation of Jo becomes 


unreliable [S6]. 

For our present problem, it is important to capture 
the gradual increase of the spin disorder with increas¬ 
ing temperature leading to decreased exchange splitting 
and band broadening. It is likely that these effects are 
not sensitive to the moderate uncertainties in the calcu¬ 
lation of Jq . The latter only affect the relative degree 
of disordering for Fe and Co, while the most interest¬ 
ing anomalies in K{x^T) occur close to pure Fe 2 B and 
C 02 B. Therefore, we adopt the following scheme based on 
the values of Jq calculated above. The MFA equations 
for a two-component alloy described within the Heisen¬ 
berg model contain four component-resolved parameters 
J^jy defined as the exchange interaction of the spin of 
atom type /i with atoms of type v everywhere else in the 
crystal. Introducing the pair exchange parameters 
we have and 

where Xy is the concentration of the component u. In 
our case this translates to J^ = ~ x)Jyf + xJfc and 

Jo =J1 — x)Jfc + xJeci where we took into account 
that Jfc = (F stands for Fe and C for Co). We 
further fit the concentration dependence of Jpc and Jee 
to a linear and Jff to a quadratic polynomial in x so as 
to best approximate the concentration dependence of J^ 
and Jq in Fig. SI. This fitting gives Jff = 406 + 220x^, 
Jfc = 225 — llx, and Jee = 116 — 41x in meV units. 
(The linear term in Jff is negligibly small.) 

The coefficients in the distribution functions for 
spin orientations are determined from the solution of the 
MFA matrix equation Ta = Jm, where d is a column 
vector with two elements J is a 2 x 2 matrix with el¬ 
ements J^jy, and m is a column vector (of reduced compo¬ 
nent magnetizations) with elements x^L{a^)^ L{a) being 
the Langevin function. The Curie temperature Tc cor¬ 
responds to the vanishing of both a^. The MCA energy 
K is calculated as a function of x and TjTc^ To facili¬ 
tate the comparison with experimental data, the iF(x, T) 
curves in Fig. 1 of the main text are plotted using the ex¬ 
perimental values of Tc for each concentration [S7]. 

In the (Coo. 9 NiQ.i) 2 B system the spin moments of Ni 
atoms are quite small, and in principle they should not be 
treated within the DLM model. However, such treatment 
does not introduce significant errors, because disorder of 
these small spin moments has a negligible effect on the 
electronic structure. Therefore, in this system we treated 
both Co and Ni within the DLM model with identical 
values of a for Co and Ni. In this way the system auto¬ 
matically tends to a paramagnetic state at a ^ 0, and 
the self-consistent calculation of the local moment of Ni 
at each temperature is avoided. 
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SPIN-ORBIT COUPLING AND 
MAGNETOCRYSTALLINE ANISOTROPY 

The spin-orbit coupling (SOC) is included as a per¬ 
turbation of the potential parameters in the Green’s 
function-based tight-binding linear mufhn-tin orbital 
(LMTO) method [S8, S9], and the resulting band prob¬ 
lem is solved exactly. This is analogous to the so-called 
“pseudoperturbative” treatment of SOC in conventional 
band-structure methods, in which the eigenvalues of the 
Hamiltonian perturbed by SOC are found exactly [SIO- 
S12]. The MCA energy is calculated as the difference in 
the single-particle energy for two directions of the mag¬ 
netization, M II X and M || z, calculated with the same 
potential parameters and distribution functions for spin 
orientations. 

The coherent-potential equations involve an integra¬ 
tion over the orientations of the spin local moment on 
each magnetic atom treated within the DLM model. For 
collinear magnetic orderings in the non-relativistic case 
the axial spin symmetry with respect to the direction of 
the magnetic order parameter is retained, and the in¬ 
tegration over the azimuthal angle can be handled an¬ 
alytically [SI]. SOC breaks this symmetry, and a full 
integration over the sphere needs to be taken. Here we 
used an 88-point quadrature for this integration, which 
was found to provide well-converged results. Thus, for¬ 
mally we apply the coherent potential approximation to 
an 176-component alloy (88 orientations of the local mo¬ 
ment for Fe and Co). Each such component is described 
by LMTO potential parameter matrices which are first 
calculated in the reference frame of the local moment 
and then rotated to the prescribed direction by a rota¬ 
tion operator generated by the total angular momentum 
operator J [SI]. 

In the analysis of the underlying mechanisms of MCA, 
we employ two approximate decompositions of K. The 
first one utilizes the calculated anisotropy of the SOC 
energy: 


The second decomposition is defined as follows. First, 
we can write the single-particle energy, calculated with 
SOC included, for magnetization direction n as 

TTin TTin 

ti/ P th p 

E%= J ENn{E)dE = E°pQ,ai + j{E - E%)N^{E)dE 

( 2 ) 

where Qvai is the total valence charge, N^{E) the density 
of states (DOS), and E'p can be set to the value of the 
Fermi energy calculated without SOC or to Ep for some 
specific orientation of n. The replacement of E^ by E^ 
in the upper limit of the last term introduces an error 
5E^ - N{Ef){E^ - E%f. We now define 

K, = J{E- E%)AN,{E)dE, (3) 

where Ncr{E) = — ^lmTrG(E) is the DOS of spin a. 
The sum -f differs from K only in the term 6Ex — 
6Ez^ which we have verified to be negligible. 

The definition (3) is used explicitly in the calculations 
of Ka. It is, however, useful to observe that in the per¬ 
turbative regime there is a relation between and ATcra' 
at T = 0. This can be seen by following the derivation in 
Ref. S8 while sorting out the spin-dependent terms. The 
second-order correction to the density of states is 

SN,iE) = --lmTvGoaVGoVGoa (4) 

TT 

where V is the perturbing SOC operator and Go the spin- 
diagonal Green’s function calculated without SOC. Using 
the cyclic property of the trace and the relation Gq^ ~ 
—dGoaldE^ which in CPA is satisfied approximately as 
long as disorder is not too strong [S8], we obtain 

= (5) 

a' 

where 


= -^A^Im / V,,,G,,,dE 


(E - E%) 


-VGoa’VdE. 


Here and in the following we use the notation AA to 
denote the difference in the values of a quantity A for 
magnetization oriented along x and along z. U is the 
perturbing SOC operator, and G the Green’s function 
calculated with SOC. The quantity Kso is naturally sep¬ 
arated in four spin contributions K^a' using the identity 
2(SL) = - {L,,)u + {L+)u + {L-)n [S8]. The 

analysis of Kso is useful, because it approximates K 
well unless second-order perturbation theory is strongly 
violated [S13]; in (Fei_a,Coa,) 2 B the concentration de¬ 
pendence of Kso agrees very well with that of K. This 
decomposition, however, loses its utility at finite temper¬ 
atures. 


On the other hand, inserting the first-order correction for 
Gcrcr' in (1), we have 

K^a' ~ j V^^>Goa'V^>^GoadE. (7) 

Note that, while it can be employed in a perturbative 
calculation of MCA energy [S12], Eq. (7) is not used ex¬ 
plicitly here. Using integration by parts and ignoring the 
relatively small energy dependence of the SOC constants, 
it is now easy to show 




Combining this with (5), we find, in the perturbative 
regime, 

(9) 

where a ^ a. The relations (8) hold when the pertur¬ 
bative approximations are admissible, i.e., as long as the 
quasi-degenerate states near do not dominate in the 
MCA energy [S12, S14]. 

If spin-off-diagonal band mixing, coming from the 
L±S^ terms in the SOC operator, can be neglected, then 
K(j ^ Ka-a comes entirely from the mixing of states in 
spin channel a by the L^Sz operator, and K^j^- represents 
a perturbative approximation for Kcr- This correspon¬ 
dence holds at all concentrations where is small, i.e., 
sufficiently far from x = 1 [S8]. However, even when 
is appreciable, the reciprocal-space resolution of K^j can 
be used to identify the hot spots in reciprocal space and 
their origin. Indeed, mixing of nearly-degenerate bands 
of the same spin a near Ef generates hot spots in 
while strong off-diagonal band mixing should give rise 
to hot spots in both and K^. The analysis of Kcr 
naturally extends to finite temperatures. 

COMPARISON OF K{T) WITH EXPERIMENT 

The calculated K(T) curves for Fe 2 B and C 02 B are 
compared with experimental data in Fig. S2. The ex¬ 
perimental temperature dependence is well reproduced 
in the calculations. The differences between different ex¬ 
periments are discussed in Ref. S16. 

SPECTRAL FUNCTIONS FOR FesB 

Fig. S3 shows the Bloch spectral functions in Fe 2 B at 
T = 0 and TjTc = 0.7. These figures can be com¬ 
pared with Fig. 2 of the main text corresponding to the 
(Feo.o 5 Coo. 95 ) 2 B composition. At the Co-rich end the 
dominant effect of spin disorder is the Stoner-like re¬ 
duction of exchange splitting with only moderate band 
broadening. In contrast, at the Fe-rich end there is very 
strong band broadening, particularly in the 1 eV energy 
window below the Fermi level. 


MAGNETOCRYSTALLINE ANISOTROPY IN 
REAL SPACE 

The main text presents the analysis of MCA in re¬ 
ciprocal space. In principle, a real-space analysis could 
provide an alternative description. Fig. S4 compares the 
term with the estimated single-site contribution to 
it. The single-site terms for Fe (or Co) were computed by 
setting the SOC parameters to zero for all atoms except 




FIG. S2. Comparison of theoretical and experimental K(T) 
curves for (a) Fe2B, (b) C02B. Solid black lines: theory; red 
dotted lines: Ref. S15; green dashed lines: Ref. S16. 


Fe (or Co) on one of the four transition-metal sites in the 
unit cell. If MCA were dominated by single-site terms, 
the concentration-weighted average of these terms would 
coincide with but Fig. S4 shows not even a corre¬ 
lation between them. Clearly, reciprocal space analysis 
is preferable to the real-space decomposition of MCA in 
this itinerant system. 
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